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Abstract. We have derived an empirical calibration of the abundance of S/H as a function of the S23 parameter, defined using 
the bright sulphur lines of [SII] and [SIII]. Contrary to what is the case for the widely used O23 parameter, the calibration 
remains single valued up to the abundance values observed in the disk HII regions. The calibration is based on a large sample 
of nebulae for which direct determinations of electron temperatures exist and the sulphur chemical abundances can be directly 
derived. ICFs, as derived from the [SIV] 10.52 yu emission line (ISO observations), are shown to be well reproduced by Barker's 
formula for a value of or = 2.5. At any rate, only about 30 % of the objects in the sample require ICFs larger than 1.2. The use 
of the proposed calibration opens the possibility of performing abundance analysis with red to IR spectroscopic data using S/H 
as a metallicity tracer. 
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1. Introduction 

Oxygen is the main abundance tracer in HII regions and HII 
galaxies, but its abundance is rather uncertain in those cases in 
which no direct determinations of the electron gas temperature 
exist therefore requiring the use of empirical or semi-empirical 
methods. These methods are based on the cooling properties of 
ionised nebulae which ultimately reflect on a relation between 
emission line intensities and oxygen abundance. In fact, when 
the cooling is dominated by oxygen, the electron temperature 
depends inversely on oxygen abundance. Since the intensities 
of coUisionally excited lines depend exponentially on tempera- 
ture, a relation is expected to exist between these intensities and 
oxygen abundances. The O23 parameter, also known as R23 and 
defined as the sum of the intensities of the [Oil] AA 3727,29 A 
and [OIII] AA 4959, 5007 A emission lines relative to U/3 (Pagel 
et al., 1979), has been widely used for these purposes. 

The relation between O23 and oxygen abundance is how- 
ever two-folded since at high metallicities, the efficiency of the 
oxygen as a coohng agent decreases the strength of the oxygen 
emission lines while at low metallicities the cooling is mainly 
exerted by hydrogen and the oxygen line strengths increase 
with metallicity. A value of 12+log(0/H) of about 8.2 divides 
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the two different abundance regimes. Although some line ratios 
have been proposed in order to break this degeneracy, the fact 
that a large number of HII regions/HII galaxies lie right on the 
turnover region is of great concern. 

On the other hand, the use of S/H as an abundance tracer 
has been frequently overlooked. Similarly to oxygen, sulphur 
is an element produced in massive stars through explosive nu- 
cleosynthesis and its yield should follow that of O. Nebular S/H 
abundances are therefore expected to follow O/H and the S/O is 
expected to remain constant at approximately the solar neigh- 
borhood value, log S/O ^ -1.6 (e.g. Lodders 2003; Bresolin et 
al. 2004). Empirical tests exploring this have been performed 
confirming the cosmic nucleosynthetic ratio though the results 
of some works suggest that this relation should be explored 
further, particularly at the not well known metallicity ends: ex- 
tremely metal deficient HII galaxies (i.e. very low O/H) and 
HII regions in the inner disk of galaxies (i.e. metal rich central 
parts with highest O/H abundances). 

The strong nebular lines of sulphur are analogous to those 
of oxygen and hence similar reasonings may be put forward re- 
garding the use of the S23 parameter, S23 = [SII]/l/l6717,6731 
+ [SIII]/l/l9069,9532)/Hy8(Vi'Ichez & Esteban, 1996) as an al- 
ternative abundance indicator (Diaz & Perez-Montero, 2000), 
which in fact presents several advantages against oxygen: 1) 
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Table 1. Reddening corrected line fluxes for sulphur normalized to I(Hy6)=100 measured in spectra of SDSS low metallicity 
galaxies referred in Kniazev et al. (2003), along with their reddening constants 



Object [5///]^6312A [SII]A6116A [SII]A6131A [SIII]A9069A C{H/3) 



SDSSJOm+ 1342 


0.9 = 


tO.2 


9.9 ± 


0.2 


8.0 H 


tO.4 


5.8 ± 


0.4 


0.03 


KU G0203 - 100 






23.0 H 


= 1.1 


18.6 


± 1.1 


19.6=1 


: 1.1 


0.19 


//50822 + 3542 


1.4 = 


tO.3 


2.8 ± 


0.2 


2.2 ± 0.2 


6.0 ± 


0.2 


0.17 


/Zwl8-A'W 


0.8 = 


tO.l 


2.0 ± 


0.3 


1.4 ±0.1 


3.9 ± 


0.3 


0.05 


/Zwl8-5£ 


0.5 = 


tO.2 


3.9 ± 


0.2 


2.8 H 


tO.2 


3.4 ± 


0.5 


0.12 


5 S5 1102 + 606 


1.6 = 


tO.3 


9.6 ± 


0.6 


6.7 H 


tO.6 


10.1 d 


:0.6 


0.15 


A1116 + 517 


1.3 = 


tO.3 


7.9 ± 


0.2 


5.3 H 


tO.2 


7.7 ± 


1.8 


0.17 


5D5'5'yil21 + 0324 


1.3 = 


tO.4 


12.2 H 


= 0.8 


9.6 H 


tO.8 


4.2 ± 


0.6 


0.20 


505571201 + 0211 


1.0 = 


tO.2 


5.4 ± 


0.2 


3.7 ± 0.2 


5.0 ± 


0.3 


0.14 


CGCG269 - 049 


1.4 = 


tO.3 


6.5 ± 


0.3 


3.1+0.3 


5.9 ± 


0.7 


0.10 



due to the longer wavelengths of the lines implied, its relevance 
as a cooling agent starts at lower temperatures (higher metal- 
licities) what makes the relation to remain single-valued up to 
solar abundances; 2) their lower dependence on electron tem- 
perature, renders the lines observable even at over-solar abun- 
dances; 3) the lines of both [SII] and [SIII] can be measured 
relative to nearby hydrogen recombination lines, thus minimiz- 
ing the effect of any reddening and/or calibration uncertainties. 
On the negative side, [SIII] lines shift out of the far red spec- 
tral region for redshifts higher than 0. 1 . Given these properties, 
here we reccommend the use of spectroscopy in the red-to-near 
infrared wavelength range in order to derive physical properties 
and abundances, of HII regions, including the lines from [SIII]/l 
6312 A, [SII]/16717,31 Auptoflie [SIII]^9069,9532 A. These 
lines can constitute a convenient analogue in this wavelength 
range to the role played by the [Oil] and [OIII] lines in the 
optical. In additon, within this wavelength range the effects of 
the extinction are less severe, thus favouring the application of 
this calibration to the study of HII regions located in the Milky 
Way disk, in the inner regions of galaxies, or to star-forming 
galaxies/regions affected by a large amount of extinction. 

In this paper, following earlier work by Christensen et al. 
(1997) and by Vermeij et al. (2002), we have explored the be- 
haviour of the sulphur emission lines in order to provide a use- 
ful abundance calibration of S/H versus S23 over the whole 
abundance range. This empirical calibration is based on the 
bright sulphur lines and encompasses the whole range of abun- 
dance currently found in HII region studies. This calibration 
appears to be independent of both, 0/H and O23, and is espe- 
cially suited for observations including (only) the red to near- 
infrared spectral range (from [SIII] /I 6312 A, Ho- up to the 1/v 
CCD cut-off). For this wavelength regime, the sulphur lines can 
be easily scaled to their respective closest Balmer or Paschen 
recombination line, H/l -instead of to Wfi- and S 23 can be writ- 
ten as follows: 



where S'23 is the scaled value of S23 and (HXjHP)o is the the- 
oretical case B recombination ratio. 



The calibrations can be extended further to the infrared to 
include the [SIV] 10.52 pL line to S23 and defining a new param- 
eter, S234 (Oey & Shields, 2000). However, the lack of [SIV] 
data renders difficult to check its reliability. 

In the next section we present the sample of objects that we 
have used here to stablish the sulphur abundance calibration as 
well as all the selected objects for which reliable ISO [SIV] 
data have been published. Section 3 is devoted to the study of 
the ionisation correction factor (ICF),-not yet well stablished 
for sulphur-, using selected data of HII regions including [SIV] 
infrared line fluxes, together with predictions of photoionisa- 
tion models. In Section 4 we present and discuss the proposed 
empirical calibration of S/H versus S23 and we summarize our 
conclusions. 

2. The selected sample of objects and sulphur 
abundances 

Our sample is a combination of different emission line ob- 
jects ionised by young massive stars: diffuse HII regions in 
the Galaxy and the Local Group, Giant Extragalactic HII re- 
gions (GEHR) and HII galaxies, and therefore does not in- 
clude planetary nebulae or objects with non-thermal activity. 
The emission line data along with their corresponding errors 
have been taken from the literature for all objects except for the 
extremely metal-poor galaxies studied by Kniazev et al. (2003) 
that present the strong [SIII] emission line at 9069 A. The spec- 
tra of these objects have been taken from the data base of the 
DR3 of the Sloan Digital Sky Survey (SDSS') and the Une 
fluxes have been measured and analyzed using the task SPLOT 
of the software package IRAF ^, following the same procedure 
as in Perez-Montero & Diaz (2003; hereinafter PMD03). We 
have found a good agreement between the oxygen emission 
lines listed by Kniazev et al. (2003) and our measurements. The 
reddening corrected sulphur emission line fluxes, normalized to 
I(Hj8) = 100, which are relevant for this work are listed in table 
[2 together with their corresponding reddening constants and 
errors. 

' The SDSS Web site is |http://www.sdss.OTg| 
- IRAF is distributed by the National Optical Astronomy 
Observatory 
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Table 2. Sulphur abundances and values of log S23 for the sample of objects. Objects marked with "('') have S^(S^^) abundances 
derived using electron temperatures based on models from t([OIII]). See text for details. 



Ob ject 


Ref." 


12 + log (^) 


log 5 23 


Ob ject 


Ref." 


12 + log(^) 


log 5 23 






M51 - CCMIO 


1 


6.95 + 0.10° 


0.02 ± 0.05 


N4A 


4 


6.74° 


0.02 


M51 - CCM53 


1 


6.95 ± 0.14° 


01 + 06 


N19E 


4 


6.59"'* 


0.00 


M51 - rrM'54 


1 


6 91 +012" 


08 + 05 


iV191A 


4 


6.80" 


0.14 


M51 - CCM55 


1 


7.04 ±0.11" 


05 ± 06 


Af80 


4 


6.02"'* 


-0.25 


M51 - CCM51 


1 


6 82 + 15° 


-0 03 + 06 


iV83 


4 


6.24"'* 


-0.09 


M5 1 - CCMll 


1 


5 90 + 13° 


-0.03 ± 0.06 


yV13 


4 


6.20 


-0.08 


NGC1232 - 04 


2 


6.43 ± 0.21° 


-0 05 + 08 


JV81 


4 


6.29* 


-0.07 


iVGC 1232 - 05 


2 


7 01 + 16° 

/ •\_/X _!_ ±\J 


16 + 17 


Af66 


4 


6.01""* 


-0.31 


NGCMM - Cn 


2 


7 02 + 21° 


08 + 06 


Wf7r604 - A 


5 


fi 82 + n 03"'* 


00 + 02 


NGC1232 - 14 


2 


6.84 ± 0.16° 


0.16 ± 0.07 




5 


6 78 + 02" * 


09 + 03 

Kj y _!_ yj .yj.^ 


NGC 1365 - 05 


2 


6.92 + 0.08° 


0.08 ± 0.01 


NGC604 - D 


5 


7.03 + 0.13° 


0.18 ± 0.01 


iVGC1365 - 08 


2 


6 81 +0 08° 


07 + 01 


NGC6Q4 - E 


5 


7 00 + 14° 


23 + 03 

\jt£^^ _!_ yj *yj .^y 


NGC1365 - 14 


2 


6.96 ±0.11" 


0.05 ± 0.01 


NGC604 


5 


6 85 + 01"'* 


06 + 01 


AfGC1365 - 15 


2 


7 00 + 05° 


23 + 01 


NGC604 


8,5 


5 79 + 09* 


-0 01 + 07 

yJ ^yj \. _!_ yj ^yj 1 


NGC1365 - 16 


2 


6.69 ± 0.13° 


03 + 01 


NGC604 


8, 22 


6 66 + 16 


-0 01 + 07 


NGC2997 - 04 


2 


6.85 ± 0.18° 


18 + 08 

V/.Xvj _!_ \J u\Jyj 


NGC604 


8, 33 


6 68 + 09* 

yj^yjij _!_ yj ,yj y 


-0 01 +0 07 


NGC29'il - 05 


2 


fj 07 + 1 6° 


1 8 + 08 


Af<7r595 


8, 33 


6 78 + 15* 

\Jt 1 \j -i- \j, i..y 


13 + 08 


NGC2991 - 06 


2 


6 92 + 1 6° 


14 + 07 


/C131 


8, 22 


6 50 + 1 8* 


06 + 08 


NGC2991 - 07 


2 


6.92 ± 0.17" 


13 + 08 


/C131 


8, 33 


6.65 +0.15* 


0.06 ± 0.08 


NGC2991 - 13 


2 


7.10 ± 0.14" 


0.15 ± 0.07 


WGC588 


8, 22 


6.42 + 0. 14* 


0.08 ± 0.07 


NGC5236 - 03 


2 


6.92 ± 0.17° 


0.12 ± 0.07 


WGC588 


8, 33 


6 63 + 09* 

—1— yj .yj y 


08 + 07 

yj^yjyj _!_ yjtyji 


A?(7r5236 - 06 


2 


7 00 + 20° 


1 2 + 08 

yjt xzj _!_ 


NGC 541] 


8,24 


6 18 + 13* 


-0 1 7 + 07 

yjt X / _!_ yj ^yj 1 


NGC5236 - 11 


2 


7.35 ± 0.1 r 


30 + 07 


/CIO - 2 


8,23 


6 61 + 06* 

_L_ V/ - V/ y.t 


09 + 06 


A/f7r5236 - 16 


2 


7.25 ± 0.13° 


25 + 07 


//7w40 


8, 7 


n + 05 

\j t L .J _!_ yj ,y}^ 


_0 22 + 06 


NGC62& - H13 


3 


6 41 + 07° 


-0 05 + 02 


//7w40 


8, 14 


6 06 + 12 

yj , _!_ V/ . X ^ 


-0 22 + 06 


NGC1232 - CDTl 


3 


7.13 ± 0.16" 


11 +0 03 

\/« X X -i- v/ ■ v/».7 


nzw40 


8, 21 


6 11 +0 06* 

\i» X X _L_ yj • \jyj 


-0 22 + 06 


NGCM31 - CDTl 


3 


6 52 + 1 6" 


21 + 03 


TI7,w40 

X X ' ' YV^xJ 


8,23 


6 14 + 06* 

Vi.XT^ \J,\J\J 


-0 22 + 06 


NGC1232 - CDT3 


3 


6 92 + 09° 

\j*y^ _!_ \j*\jy 


16 + 02 


n7w40 

X 1 ' ■ r K v/ 


8, 24 


6 12 + 06* 

_!_ yj -yjyj 


-0 22 + 06 


NGC 1232 - CDT4 


3 


6.93 ± 0.07° 


0.16 + 0.02 


A'GC4861 


8,6 


6.37 + 0.12 


-0.02 ± 0.12 


NGC2467 


4 


6.39* 


-0.03 


WGC2363 


8, 26 


5 83 + 14 


-0 40 + 13 




4 


6.99 


0.19 


IZwl23 


8,7 


5 92 + 16* 


-0 27 + 15 


Mil 


4 


6.92* 


0.20 


Mrk36 


8, 15 


5.91 + 0.17* 


-0.19 ± 0.17 


M20 


4 


6.99° 


0.20 


Mrk36 


8, 7 


5.88 + 0.14* 


-0.19 ± 0.17 


NGC3576 


4 


7.14* 


0.46 


Irl f fVyj\./\J 


8, 24 


6 00 + 14* 

yj.y/Vf _!_ y/. x^ 


-0 34 + 13 


Orionl 


4 


6.93 


0.03 


I Zw 18 


8,23 


5.49 + 0.15"'* 


-0.68 ±0.18 


Orion2 


4 


6.97* 


0.26 


IZwlS 


8,24 


5.51+0.19"'* 


-0.68 ±0.18 


N59A 


4 


6.60° 


-0.01 


MlOl -//681 


9 


6.06 ± 0.09" * 


-0.20 ± 0.07 


N44B 


4 


6.84° 


0.09 


M51 - CCMIO 


10 


7.16 ±0.20° 


0.20 ± 0.05 


N55A 


4 


6.88° 


0.04 


M51 -CCM72 


10 


7.56 ± 0.02° 


0.41 ± 0.00 


NIUD 


4 


6.90" 


0.25 


AfGC2403- VS35 


11 


6.77 ±0.17 


0.10 ±0.03 


NniA 


4 


6.83 


0.12 


AfGC2403 - 7524 


11 


7.04 ±0.11 


0.10 ±0.04 


N\59A 


4 


6.53°-* 


0.02 


A^GC2403 - 1/538 


11 


6.87 ± 0.09 


0.06 ± 0.02 


N2UC 


4 


6.75° 


0.12 


iVGC2403 - VS44 


11 


6.63 ± 0.07 


0.04 ± 0.02 



References are: 1. Bresolin et al., 2004; 2. Bresolin et al., 2005; 3. Castellanos et al., 2002; 4. Dennefeld & Stasinska, 1983; 5. Di'az et al., 1987; 6. Dinnerstein & Shields, 1986; 
7. French, 1980; 8. Garnett, 1992; 9. Gamett & Kennicutt, 1994; 10. Gamett et al., 2004; 11. Gamettetal., 1997; 12. Gonzalez-Delgado et al., 1995; 13. Gonzalez-Delgado et al., 
1994 14. Guseva et al., 2000; 15. Izotov & Thuan, 1988; 16. Izotov et al., 1994; 17. Izotov et al., 1997; 18. Kennicutt et al., 2003; 19. Kinkel & Rosa, 1994; 20. Kniazev et al., 
2003; 21, Kunth & Sargent, 1983; 22. Kwitter & Aller, 1981; 23. Lequeixx et al., 1979; 24, Pagel et al., 1992; 25. Pastoriza et al., 2003; 26. Peimbert et al., 1986; 27. P&ez-Montero 
& Dfaz, 2003; 28. SkiUman & Kennicutt, 1993; 29. SkiUman et al., 1994; 30. Terlevich et al., 1991; 31. Vermeij et al., 2002; 32. Vflchez & Esteban, 1996; 33. Vflchez et al., 1988; 
34.This work. 



For all the objects of the sample measurements of the emis- 
sion lines of [SII] at 6717,31 A and of [SIII] at 9069,9532 A 
exist, thus allowing the simultaneous determination of the S23 
parameter and the abundances of S"^ and S^"^. 



The physical conditions of the ionised gas in each sam- 
ple object, including electron temperatures, electron density 
and sulphur abundances, have been computed from the origi- 
nal emission line data using the same procedures as in PMD03, 
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Table 2. Continued. 



Ob ject 


Ref." 


12 + log (^) 


log 5 23 


Ob ject 


Ref." 


12 + log (^) 


log 5 23 






NGC2403 -VS51 


11 


6.75 ±0.17 


0.06 + 0.06 


NGC3310 - M 


25 


6.85 + 0.01* 


0.06 ± 0.01 


NGC2403 -VS3 


11 


6.62 ± 0.04 


0.06 ± 0.03 


IIZw40 


27,7 


6.06 ± 0.04 


-0.28 ± 0.02 


NGC2403 - VS49 


11 


6.59 ± 0.17 


0.06 ± 0.05 


IIZw40 


27, 14 


6.04 ± 0.05 


-0.28 ± 0.02 


NGC2403 - VS48 


11 


6.60 ± 0.26 


0.08 ± 0.09 


llZw40 


27,21 


6.06 ± 0.04 


-0.28 ± 0.02 


NGCnU-A 


12 


6.81 ± 0.05" 


0.19 ±0.01 


IIZw40 


27,23 


6.10 ±0.04 


-0.28 ± 0.02 


NGCnU-NWQ 


12 


6.99 ±0.11" 


0.18 ±0.05 


IIZw40 


27,24 


6.06 ± 0.04 


-0.28 ± 0.02 


NGCllU - B 


12 


6.31 ±0.06" 


0.01 ± 0.02 


Mrk5 


27, 15 


6.19 ±0.09 


0.07 ± 0.02 


NGCll 14 - C 


12 


6.41 ±0.18" 


0.00 ± 0.09 


5BS0749 + 568 


27,17 


6.04 ±0.11 


-0.17 ±0.04 


NGCnU-N 


12 


6.95 + 0.11" 


0.22 ± 0.05 


5 S5 0926 + 606 27,17 


6.04 ±0.11 


-0.17 + 0.04 


NGC2363 - A2 


13 


5.92 ± 0.08 


-0.41 ± 0.02 


MrA:709 


27,30 


5.97+0.10 


-0.04 + 0.02 


MlOl -moi3 


18 


7.12 ±0.07 


0.21 ± 0.02 


Mrk22 


27,16 


6.12 + 0.11 


-0.11+0.05 


MlOl -HUOS 


18 


6.95 ± 0.06 


0.18 ± 0.02 


Mrkl434 


27, 17 


5.76 ± 0.08 


-0.38 ± 0.03 


MlOl -//1159 


18 


6.68 ± 0.09 


0.10 ±0.03 


Mrk36 


27,7 


5.96 ±0.12 


-0.21 ± 0.04 


MlOl -Hino 


18 


7.08 ± 0.08 


0.36 ± 0.03 


Mrk36 


27,15 


5.98 ±0.13 


-0.21 ± 0.04 


MlOl -Hine 


18 


6.82 ± 0.06 


0.14 ±0.02 


VIIZw403 


27, 17 


6.10 ± 0.11 


-0.27 ± 0.04 


MlOl -m2l6 


18 


6.55 ± 0.07 


0.01 ± 0.03 


UM461 


27, 15 


5.74 ± 0.06 


-0.35 ± 0.02 


MlOl -H\2% 


18 


6.81 ± 0.07 


0.11 ± 0.02 


UM462 


27, 15 


6.13 ± 0.10 


-0.13 ± 0.04 


MlOl -H143 


18 


6.73 ± 0.09 


0.13 ±0.03 


Mrk209 


27, 17 


5.94 ±0.10 


-0.28 ± 0.04 


MlOl - H149 


18 


6.72 ± 0.07 


0.12 ±0.02 


IZwli-SE 


28 


5.31 ±0.12* 


-0.78 ± 0.07 


MlOl -H336{S5) 


18 


7.10 ± 0.09 


0.21 ± 0.03 


IZwlS - NW 


28 


5.13 ±0.10* 


-0.93 ± 0.08 


MlOl -H4Q9 


18 


6.63 ± 0.07 


0.08 ± 0.03 


UGC44S3 


29 


5.77 ± 0.13 


-0.39 ± 0.05 


MlOl - H61 


18 


6.60 ± 0.14 


0.07 ± 0.03 


N160-A1 


31 


6.91 ± 0.06 


0.17 ±0.03 


MIQI-NGC5411A 


18 


6.29 ± 0.07 


-0.11 ±0.03 


N160-A2 


31 


6.86 ± 0.07 


0.14 ±0.03 


MlOl - NGC5411B 


18 


6.31 ±0.08 


0.02 ± 0.03 


N159 - 5 


31 


7.01 ±0.10 


0.22 ± 0.03 


MlOl -NGC541\C 


18 


6.42 ± 0.08 


-0.05 ± 0.03 


N157 -B 


31 


6.87 ± 0.09 


0.24 ± 0.04 


MlOl -7VGC5471D 


18 


6.51 ± 0.10 


-0.02 ± 0.03 


30Dor - I 


31 


6.72 ± 0.08 


0.06 ± 0.02 


MlOl -S5 


19 


7.31 ± 0.04 


0.31 ± 0.00 


30Dor - 2 


31 


6.80 ± 0.09 


0.11 ± 0.02 


SDSSJQ133+ 1342 


20, 34 


5.58 ±0.13" 


-0.42 ± 0.03 


30Dor - 3 


31 


6.87 ± 0.06 


0.14 ±0.02 


KUG()2()3 - 100 


20, 34 


6.02 ± 0.04" 


0.04 ± 0.05 


30Dor - 4 


31 


6.79 ± 0.08 


0.12 ±0.02 


//50822 + 3542 


20, 34 


5.27 ±0.19" 


-0.59 ±0.01 


Nil -A 


31 


6.78 ± 0.09 


0.10 ± 0.03 


/Zwl8 - NW 


20, 34 


5.08 ±0.18" 


-0.78 ± 0.04 


NS3-B 


31 


6.82 ± 0.07 


0.13 ± 0.03 


IZwl%-SE 


20, 34 


5.19 + 0.23" 


-0.73 ± 0.06 


N19-A 


31 


6.96 ± 0.07 


0.19 ±0.02 


5S5 1102 + 606 


20, 34 


5.67 + 0.19" 


-0.29 ± 0.02 


N4-A 


31 


6.79 + 0.07 


0.10 ± 0.02 


A1116 + 517 


20, 34 


5.49 + 0.31" 


-0.40 ± 0.09 


N%%-A 


31 


5.92 + 0.10 


-0.32 ± 0.03 


SDSSJim +0324 


20, 34 


5.49 ±0.31" 


-0.44 ± 0.06 


N66 


31 


6.23 ± 0.15 


-0.08 ± 0.05 


5D5571201 + 0211 


20, 34 


5.28 ±0.18" 


-0.58 ± 0.02 


NSl 


31 


6.22 ±0.10 


-0.17 ±0.03 


CGCG269 - 049 


20, 34 


5.39 ±0.21" 


-0.50 ± 0.05 


5209 


32 


5.76 ± 0.20* 


-0.35 ± 0.07 


MlOl - NGC5455 


22 


6.75 ±0. 10"'* 


0.14 ±0.05 


5127 


32 


6.42 ± 0.29* 


-0.12 ±0.07 


MlOl -NGC5411 


22 


6.31 ±0.10"* 


-0.06 ± 0.05 


5128 


32 


6.40 ± 0.29* 


-0.09 ± 0.06 


NGC3310-A 


25 


6.82 ± 0.03" 


0.00 ± 0.02 


MA2 


33 


6.94 ± 0.34* 


0.08 ± 0.02 


NGC3310-B 


25 


6.78 ± 0.02* 


0.09 ± 0.03 


yVGC604 


33 


6.76 ± 0.03* 


0.06 ± 0.01 


NGC3310-C 


25 


7.03+0.13 


0.18 ±0.01 


AfGC595 


33 


7.02 ± 0.09* 


0.36 ± 0.02 


NGC3310-E 


25 


7.00 + 0.14 


0.23 ± 0.03 











References are: 1. Bresolin et al., 2004; 2. Bresolin et al., 2005; 3. Castellanos et al., 2002; 4. Dennefeld & Stasinska, 1983; 5. Dfaz et al., 1987; 6. Dinnerstein & Shields, 1986; 

7. French, 1980; 8. Garnett, 1992; 9. Garnett & Kennicutt, 1994; 10. Garnett et al., 2004; 11. Garnett et al., 1997; 12. Gonzalez-Delgado etal., 1995; 13. Gonzalez-Delgado et al., 
1994 14. Guseva et al., 2000; 15. Izotov & Thuan, 1988; 16. Izotov et al., 1994; 17. Izotov et al., 1997; 18. Kennicutt et al., 2003; 19. Kinkel & Rosa, 1994; 20. Kniazev et al., 
2003; 21. Kunth & Sargent, 1983; 22. Kwitter & Alien 1981; 23. Lequeux etal., 1979; 24. Pagel et al.. 1992; 25. Pastoriza et al., 2003; 26. Peimbert et al., 1986; 27. Perez-Montero 
& Di'az, 2003; 28. Skillman & Kennicutt, 1993; 29. Skillman et al., 1994; 30. Terlevich et al., 1991; 31. Vermeij et al., 2002; 32. Vilchez & Esteban, 1996; 33. Vilchez et al., 1988; 
34. This work. 



based on the five-level statistical equilibrium model in the 
task TEMDEN and IONIC, respectively, of the software pack- 
age IRAF (De Robertis, Dufour & Hunt, 1987; Shaw & 
Dufour, 1995). The atomic coefiicients used are the same as 
in PMD03 (see Table 4 of that work). Electron densities are 



determined from the the [SII] A 6717A / A 6731A Une ratio. 
Electron temperatures have been calculated from the [SIII] (A 
9069A+^9532A)/ A 6312A Une ratio for all but 44 objects of 
the sample for which the [OUT] (A 4959A+^5007A)//l 4363A 
line ratio has been used. For these objects, marked with * in 
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table 2, a theoretical relation between [OIII] and [SIII] electron 
temperatures has been used: 

t([SIII]) = 1.05f([C»///])-0.08 

This relation is based on the grids of photo-ionisation mod- 
els described in Perez-Montero & Diaz (2005) and differs 
slightly from the empirical relation found by Garnett (1992), 
mostly due to the introduction of the new atomic coefficients 
for 8^+ from Tayal & Gupta (1999). 

Regarding [SII] temperatures, for those objects without the 
[SII] auroral lines at 4068,4074 A. we have taken the ap- 
proximation t[SII] ~ t[OII] as valid. For 124 objects of the 
sample it has been possible to derive t[OII] from the [Oil] (/I 
3726A-H/13729A) /A 7325 A Une ratio. ^ For the rest of the ob- 
jects of the sample, marked with " in table 2 , not presenting any 
auroral line in the low excitation zone, we have resorted to the 
model predicted relations between t([OII]) and t([OIII]) found 
in PMD03 that take explicitly into account the dependence of 
t([OII]) on electron density. This can affect the deduced abun- 
dances of IH^ by non-negligible factors, larger in all cases 
than the reported errors. 

For those objects for which multiple observations exist we 
have considered each one of them as independent. The final 
quoted errors in the derived quantities have been calculated by 
propagating the measurement errors in the emission lines pro- 
vided by the different authors. This information is not provided 
for the objects from Dennefeld & Stasiriska (1983; reference 4). 
The ionic abundances of sulphur, S^/H^ and S^^/H^, together 
with the values of the S23 parameter for each object are given 
in table 121 

It is important to emphasize here that this database include 
objects covering all the abundance range from low metallicity 
HII galaxies, at 1/40 of the solar abundance, up to HII regions 
in the disks of spirals populating the high metallicity range, 
up to «; 3 Zq. Although he data selected for this study have 
been obtained using different apertures and instrumental con- 
figurations and therefore do not constitute a homogeneous sam- 
ple, they have been reanalysed and ionic abundances have been 
derived in a homogeneous manner. Using this sample we can 
study on a firmer basis the empirical relationship between the 
sulphur abundance and the parameter S23. In order to do that we 
explore first the sum of the abundances of and S^^ as a func- 
tion of S23 and, on a second step, we will derive the calibration 
of the total abundance of sulphur as a function of S23. The re- 
sult is plotted in figure^ showing a relationship with very low 
scatter for which we have obtained the following quadratic fit: 

12 4- log — — — U 6.540 + 2.071 log523 +0.348 (log523)^ 

with a typical dispersion of 0.17 dex, defined as the standard 
deviation of the residuals of the points. 

3 The [Oil] /17319A+.17330A lines can have a contribution by di- 
rect recombination which increases with temperature. Using the cal- 
culated [OIII] electron temperatures, we have estimated these contri- 
butions to be less than 4 % in all cases and therefore we have not 
corrected for this effect 




-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 



log 

Fig. 1. The calculated sulphur abundance for the objects of the 
sample versus the S23 parameter The thick line represents the 
quadratic fit of the points. 

After the advent of the ISO mission, it has also been pos- 
sible to obtain data of the [SIV] at 10.52 fi emission line for 
a selected sample of objects, thus allowing the derivation of 
the S^^/H^ ionic abundance ratio. We have collected data for 
1 1 HII regions (plus one supernova remnant) in the Large and 
Small Magellanic Clouds (Vermeij et al., 2002) and for the HII 
galaxy Mrk209 (Nollenberg et al., 2002). For all these objects 
the S ionic abundances have been recalculated following the 
procedure described above and using the most recent atomic 
coefficients from Saraph & Storey (1999). The calculated abun- 
dances are listed in table 3. 

It has been a matter of debate for some time now what is the 
exact contribution of the S^^/H^ ionic abundance to the total 
abundance of sulphur. It is clear that an efficient solution to this 
problem could be reached making use of the measurements of 
the infrared lines of [SIV] -not available for a large sample of 
objects yet-. 

This question takes us straight into the issue of the ioni- 
sation correction factor, a necessary step for the derivation of 
the empirical sulphur abundance calibration. In any case, it has 
been already pointed out that a one-zone ionisation scheme 
may provide some insight into the situation for not ionisation 
bounded HII regions of moderate to low ionisation, where the 
various ions coexist throughout the nebula (e.g. Pagel 1978). 
For this reason we have presented above the general coiTelation 
between S^ + S^^ and S23. A study of the ionisation coiTection 
factor scheme for sulphur is presented in the next section. 

3. The ICF Scheme 

Perhaps one of the most difficult aspects in the derivation of to- 
tal abundances is the question of the ionisation correction fac- 
tor, ICF. The ICF for sulphur accounts for the contribution of 
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Table 3. Abundances of S^^ for the objects for which there are 
ISO observations of the [SIV] 10.52 // emission line and the 
corresponding ICF(S++S2+). 



Object Ref." 12 + log (^) ICF(S*+S^+) 



Mrk209 

Mrk209'' 
Nm-Al 
N160-A2 

A'159-5 
N157-B 
30Dor - 1 
30Dor - 2 
30Dor - 3 
30Dor - 4 

NS3-B 
N4-A 

ms-A 

N66 

mi 



5.93: 
5.68: 
6.08: 
5.90: 
6.03: 
5.54: 
6.16: 
6.24: 
6.16: 
6.15: 
5.45: 
6.17: 
5.80: 
5.73: 
5.45: 



0.20 
0.23 
0.11 
0.10 
0.10 
0.12 
0.11 
0.11 
0.10 
0.10 
0.12 
0.15 
0.10 
0.15 
0.11 



1.97 : 

1.56: 
1.15: 
1.11 : 
1.10: 
1.05: 
1.27: 
1.27 : 
1.20: 
1.23 : 
1.04: 
1.24: 
1.75: 
1.32: 
1.17 : 



0.62 
0.76 
0.05 
0.03 
0.04 
0.02 
0.10 
0.10 
0.06 
0.07 
0.02 
0.11 
0.27 
0.19 
0.07 



The conesponding references for the emission hne data of [SIV] at 10.5 fi are: 1. 
Nollenberg et al., 2002; 2. Vermeij et al., 2002 
^ Assumming a constant 5^^/.?"^ ratio from the [SIII] lines in the near-IR from 
PMD03 



the ionic species not detected in the optical. In high excitation 
nebulae a large fraction of the sulphur can be found in the S^^ 
stage, whose prominent emission lines of [SIV] are observed 
in the mid-IR at 10.52 /v. Therefore, in order to derive the total 
abundance of sulphur it is necessary to correct for the presence 
of S^"^ in those objects for which there are not observations in 
the mid-IR, as follows. 



ms) 



= ICF(S^ +S^^) 



NiH) 
_ N(S^ +S^^ +S^^) 



N{H+) 



The first proposed ICF scheme for sulphur (Peimbert & 
Costero, 1969) was based on the similarity of the ionisation 
potentials of 0+ (35.1 eV) and S^^ (34.8 eV) 



ICF(S' 



N{0) 



Nevertheless, some authors have pointed out that this rela- 
tion has a strong correlation with the ionisation degree of the 
nebula (e.g. Barker 1978; Pagel 1978), thus implying an over- 
estimation of the sulphur abundance in nebulae with low elec- 
tron temperature. Barker (1980) proposed a new relation, based 
on the photoionisation models of Stasinska (1978): 



ICF(S^ +S^^) = 



1-1- 



mo) 



-I /a 



for which he proposed a value of a = 3. The ICF from Peimbert 
& Costero corresponds to a value for a = 1 in this expresion. 



ft. 
u 



2.8 
2.6 
2.4 
2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 



Mrk 209 
Mrk 209 rev. 
O LMC 
<$> 30 Dorados 
SMC 




I a=,3 



0.0 



0.2 



0.4 0.6 
o*/o 



1.0 



Fig. 2. Ionisation correction factor for the sum of and S "^^ 
versus the 10 ratio for the sample of objects listed in table 
3. The solid lines represent the Barker's formula for different 
values of the a parameter. Squares represent different values for 
Mrk209 (see text for details), diamonds for HII regions in the 
Large Magellanic Cloud and circles in the Small Magellanic 
Cloud. 



Later, Izotov et al. (1994), based on photoionisation models 
from Stasinska (1990) gave a fit for this ICF that, in fact, is 
quite similar to the Barker formula for a -2. 

In figure 2 we show our computation of the ICF(S^ + S^^) 
as a function of log (0^/0) using the data listed in table 3, 
for three different values of a =1,2,3. It is apparent in this 
plot that most points are better matched for a values between 
2 and 3. Only one point, corresponding to a HII galaxy (Mrk 
209, solid square), shows an ICF for a value of a even lower 
than 2. Since the abundances obtained from the [SIII] line 
at 18.71 yu from the ISO observations from Nollenberg et al. 
(2002) are much higher than those obtained from the near-IR 
[SIII] (PMD03), we have recalculated S abundances assum- 
ing a constant S ^'^/S ratio. This value is showed in table 3 and 
is represented as an open square in figures 2 and 3. Although 
showing a large error bar, the new value lies within the zone of 
a between 2 and 3, in better agreement to the values predicted 
by photo-ionisation models (PMD03). 

Since the relevant lines for sulphur abundance determina- 
tion are in the red to near infrared range, it is worth trying 
to relate the ICF scheme of sulphur to the ionisation structure 
as seen in this wavelength range. This aproach has the extra 
bonus of reducing the effect of reddening in the line ratios cur- 
rently used. In order to derive a new ICF scheme for sulphur 
based only on the red-to-near infrared information, we need 
to combine available data and compute new models for ref- 
erence objects, for which we know all the relevant ionic abun- 
dances of sulphur. We have computed this relation and is shown 
in figure |3] In this plot, we present the ICF(S^ + S^^) vs. 



2.2 



2.0 



1.8 



1.6 

+ 

is 
y 

" 1.4 



1.2 



T 



Perez-Montero et al.: An empirical calibration of sulphur abundance in ionised gaseous nebulae 

8.0 



T 



T 



■ Mrk209 

□ Mrk209 rev. 

O LMC 

<i> 30 Dorados 

• SMC 

X HIIG models 



/ 

/ 

/ 

/ / 



T 



T 



T 



/50 kK 



'/ 



)kK 



4 



_L 



_L 



_L 



_L 



35 kK 
I.I.I 



_L 



_L 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 
log ([S11I]/[S11]) 



1.0 1.2 



Fig. 3. lonisation correction factor for the sum of 5 ^ and S 
versus log ([SIII]/[SII]) for the sample of objects listed in table 
3. The dashed lines represent the results from photo-ionisation 
models for CoStar single star atmospheres with the labelled 
effective temperature. Crosses represent results from photo- 
ionisation tailored models of HII galaxies. The rest of symbols 
are the same as in figure 2. 



log ([SIII]/[SII]) predicted by photo-ionisation models using 
CoStar model atmospheres of different effective temperatures 
(35kK, 40kK, 50kK, Perez-Montero & Diaz, 2005) and HII 
galaxy models (Perez-Montero & Diaz, 2004), together with 
all the observed points available. Clearly, the ICF predicted 
by a model with a single-star of Teff=35kK is always giving 
ICF=1 no matter the excitation. The ICF predicted for objects 
presenting log([SIII]/[SII]) <0.2 remains small, 1.0<ICF<1.05. 
Above log([SIII]/[SII])=0.4, the ICF model predictions begin 
to diverge. Up to some point, this behaviour seems to be fol- 
lowed by the data. Giant HII regions points cluster around the 
locus of 40kK models, except in the notable case of N88A, a 
SMC reddened, very young HII region breaking out its natal 
cloud (Heydari-Malayeri et al., 1999). In this plot this giant 
HII region is located closed to Mrk209, the only HII galaxy in 
table 3. This fact indicates a hotter ionising source and a larger 
ICF. It is suggested here that this behaviour could be consistent 
with the prediction from single burst evolutionary models (e.g. 
Stasinska et al. 2001) since typically HII galaxies host younger 
ionizing clusters than giant HII regions (Terlevich et al. 2004). 
Under this assumption, the equivalent width of Ha, an age in- 
dicator, should provide a useful constraint. 



4. Discussion 

Earlier works by Christensen et al. (1997) and by Vermeij et al. 
(2002) have explored the abundance calibration of S/H versus 
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Fig. 4. Calculated total sulphur abundance for the objects of the 
sample versus the S23 parameter The thick line represents the 
empirical cahbration proposed in this work 

S23. Christensen et al. (1997) proposed a linear S23 calibration 
as follows 

<S 



12 -Hlog(-j = 6.485 + 1.218 log523 



using data from giant HII regions and HII galaxies available to 
this date, and complementary model predictions from Stasinska 
(1990). Though these points included some high metallicity 
HII regions of M51 from Diaz et al. (1991), as well as lower 
metallicity objects from Garnett (1992), still the whole abun- 
dance range was not sufficiently well sampled. They claimed 
that more data were needed before a definite conclusion could 
be drawn. 

Vermeij et al.(2002) derived sulphur abundances for a new 
data set of optical and infrared spectra of HII regions in the 
Large and Small Magellanic Clouds. This information allowed 
these authors to derive all the ionic fractions of sulphur; includ- 
ing S^^* from ISO [SIV] 10.5 yU line for their sample objects. 
They present a comparison of their S/H abundances points vs. 
S23 and S234 with model calculations for log U = 0, -1, -2, 
and -3. The following relation was found between log S23 and 
log(S/H) 



-5.65 H- 1.50 log 5 234 



However no relation was proposed for S/H and S23. 

In this work, we present an empirical abundance calibra- 
tion, based on the bright sulphur lines, which encompasses the 
whole range of abundance currently found in HII region stud- 
ies. This calibration is firmly based on an extended homoge- 
neous data base. In figure 4 we present the relation for all the 
objects of our sample between the S23 parameter and the to- 
tal abundance of sulphur, taking into account the ionic abun- 
dances of S^, S^^ and the ICF corresponding to the formula of 
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7,0 7,5 8,0 8,5 9,0 

12 + log (0/H) 



Fig. 5. Representation of the quotient of S/O versus O/H for the 
compiled sample and the solar value 

Barker for a - 2.5 which is the value that better fits the avail- 
able points. The best quadratic fit to the data gives: 

12 + log (^j = 6.622 + 1.860 log523 +0.382 (log523f 

with a dispersion of 0.185 dex in the range of -1.0 < log(S23) < 
0.5. 

This calibration is, to first order, independent of both, O/H 
and O23, and is especially suited for observations including 
only the red to near-infrared spectral ranges (from [SIII]/l 6312 
A, Ho- up to the Ifi CCD cut-off). Though the S23 parameter is 
possibly double valued, like R23, the turnover region for S23 is 
located above the range of sulphur abundance currently found 
in disk HII regions. 

The calibration is, to some extent, affected by the ICF cal- 
culation scheme. However only 34% of the objects in our sam- 
ple show an ICF larger than 1.2 as derived from Barker's ex- 
pression (a = 2.5) and all of them show log ([SIII]/[SII]) > 0.4, 
where model predictions are more uncertain. In fact, the cali- 
bration of (S^+S^^)/H^ vs. S23 does not differ much from that 
of S/H vs. S23 (see figures 1 and 4) and the values obtained 
from both for a given S23 are within the quoted uncertainty in 
most cases. 

This calibration makes possible the use of S/H as metallic- 
ity tracer in ionized nebulae. Its translation however to an O/H 
abundance relies on the assumption that the S/O ratio remains 
constant at all abundances. This point though remains contro- 
versial (see e.g.Lodders, 2003; Bresolin et al. 2004). Figure 5 
shows the S/O ratio vs. 12+log(0/H) for the objects in our sam- 
ple. For all of them, the O/H abundance has been calculated fol- 
lowing the scheme presented in Perez-Montero & Diaz (2005). 
It can be seen from the figure that HII galaxy data are consistent 
with a constant S/O ratio, but significantly lower than the solar 
ratio. Four galaxies deviate from this trend (UGC4483, KUG 



0203-100, Mrk709 and SDSS Jl 121+0324). One of them, Mrk 
709, also shows a large value of N/O (PMD03). Regarding disk 
HII regions the dispersion is much larger and the assumptions 
of a constant S/O is highly questionable. 

Summarizing, following earlier work by Christensen et al 
and Vermeij et al, we have derived an empirical calibration of 
the abundance of S/H as a function of the S23 parameter, de- 
fined using bright sulphur lines, which we recommend as a 
useful tool to derive S/H within a large abundance range of 2 
dex, keeping a statistical error of 0.18 dexrms. This abundance 
range appears well suited to deal with objects from low metal- 
licity HII galaxies to high metallicity HII regions located in the 
inner parts of the disks of spirals. 
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